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Abstract. We study the dynamical low temperature behaviour of the Ising spin glass 
on the Bcthe lattice. Starting from Glauber dynamics we propose a cavity like Ansatz 
that allows for the treatment of the slow (low temperature) part of dynamics. Assuming 
a continuous phase transitions and ultrametricity with respect to long time scales we 
approach the problem perturbatively near the critical temperature. The theory is 
formulated in terms of correlation-response-functions of arbitrary order. They can, 
however, be broken down completely to products of pair functions depending on two 
time arguments only. For binary couplings J ~ ±1 a spin glass solution is found 
which approaches the corresponding solution for the SK-model in the limit of high 
connectivity. For more general distributions P{ J) no stable or marginal solution of 
this type appears to exist. The nature of the low temperature phase in this more 
general case is unclear. 
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1. Introduction 



Most of the work on spin glasses and other systems with frozen-in disorder is based on 
the evaluation of the free energy or the ground state energy. This comprises replica 
theory, the cavity method or the TAP-equations [1]. As an alternative stochastic 
dynamics [2] has been employed for systems with continuous freezing transition [3, 4], 
e.g. the Sherrington-Kirkpatrick (SK) model [5], and also for systems with discontinuous 
transition, e.g. the spherical spin glass with p-spin interactions [6, 7, 8]. The essential 
difference between the two approaches shows up in the thermodynamic limit. The 
computation of the free energy does not rely on any kind of dynamics and the question 
of how equilibrium states can be reached starting from non-equilibrium initial conditions 
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is not adressed. In fact diverging barriers build up in the thermodynamic hmit N ^ oo. 
As a consequence the approach via dynamics in the non-ergodic low temperature phase 
requires regularization by some long time scale, for instance the waiting time t^ after 
a quench from high temperature, by introducing slowly varying bonds [4] or slow 
cooling [9]. Typically the thermodynamic and long time limit are taken in the order 
limt^_»oo liniAT^oo- This means that diverging barriers can not be overcome or the system 
might be stuck in metastable states [8, 10]. The importance of the order in which the 
two limites are performed has been demonstrated in a different context (learning in a 
perceptron with binary synapses) [11]. For systems with discontinuous transition the 
dynamic freezing transition is higher than the temperature where replica symmetry 
breaking sets in, whereas both temperatures are identical for continuous transitions. In 
this latter case the coefficients of an expansion of the energy or the Edwards-Anderson 
order parameter near Tc are identical at least up to fourth order [9]. At zero temperature, 
however, the energy obtained within dynamics is expected to be higher than the ground 
state energy determined within replica theory. Furthermore the energy found within 
dynamics may depend on details of the cooling schedule [12]. 

Spin glasses on diluted graphs, for instance on the Bethe lattice, share finite 
connectivity with spin glasses in finite dimensions. Nevertheless they are of mean field 
character as the fully connected SK-model [5, 1]. Following earlier attempts [14, 15], the 
Ising spin glass on the Bethe lattice has been solved by Mezard and Parisi [13] using the 
cavity method on the level of one step replica symmetry breaking. The cavity method 
benefits from the local tree like structure of the Bethe lattice and deals on the level 
of IRSB with distributions of local fields and distributions of those distributions. The 
resulting functional equations have been solved numerically with a population dynamics 
algorithm for binary couplings J = ±/. 

Investigating the dynamics of systems of Ising spins one has to specify not only 
the Hamiltonian (energy) but also the kind of dynamics. In the present approach 
Glauber dynamics, i.e. stochastic single spin fiip dynamics with transition probabilities 
depending on temperature, is used. In the non-ergodic low temperature phase a 
separation of short and long time scales is assumed. The longest time scale is realized 
by a waiting time or by some other means [4, 9]. On the short time scale the 
validity of fiuctuation-dissipation theorems are assumed to hold. They are characteristic 
for equilibrium and in this context they describe equilibrium within a single valley of 
the energy landscape. On the long time scale in analogy to the cavity method [13] a 
functional equation for the distribution of histories of local fields is derived, investigating 
the iterative assembly of local subtrees of the Bethe lattice to new subtrees. 

Hierarchically connected equations are derived from this functional for correlation- 
response-functions of, in principle, arbitrary order. In order to truncate this system 
of equations, an expansion in powers of the deviation of the temperature from the 
freezing temperature is employed. Such an expansion is known for the dynamics of 
the SK-model. It turns out that correlation-response-functions of higher order can be 
broken down completely to expressions involving pair functions depending on two time 
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arguments only. For binary couplings J = ±/ a marginal solution resembling the one 
for the SK- model is found [3, 4]. 

The situation for more general distributions P{J) is unclear. There exists a critical 
temperature where the ergodic high temperature solution becomes unstable. The above 
solution does, however, not apply and no other stable or marginal solution could be 
found. The expansion contains contributions corresponding to the one step replica 
symmetry breaking solution. The original static IRSB solution [13] has been evaluated 
for J = ±/ only. It would be interesting to see whether this solution persists for general 
P(J). 

The paper is organized as follows. Sec. 2 contains a brief survey of Glauber dynamics 
and fixes some notations. In Sec. 3 the dynamics of the Ising spin glass on the Bethe 
lattice is formulated, an effective single site evolution is introduced and equilibrium 
properties are investigated. Dynamics on long time scales is dealt with in Sec. 4 including 
a formulation in terms of distributions of slow components of local fields. In view of an 
expansion valid in the neighborhood of the freezing temperature, this is reformulated 
in terms of correlation-response functions. Sec. 5 contains the expansion around the 
critical temperature in leading order for general P{J) and in next to leading order for 
J = ±J. The paper concludes with a stability analysis in Sec. 6 and a brief discussion. 



2. Glauber dynamics of a single spin 



The present investigation is based on Glauber dynamics for Ising spins. We start with 
a brief outline of the formulation [11] describing the dynamics of a single spin. The 
Hamiltonian of the spin in an external field h is 



H = -ha,. 

The actual state of the spin at time t is described by a 2-component vector 



p{t 



(V^it) \ ^ ( i(l + m(t)) 
\P-it) J [ 
Introducing the Pauli matrix 

the expectation value of the spin is given by 



1 1 



1 
-1 



(2.1) 



:2.2) 



(2.3) 



(l|a,|p(t)) =m{t). (2.4) 



It is convenient to use a notation resembling quantum mechanics, although Glauber 
dynamics is purely classical. The time evolution is described by a Liouville operator, a 
two by two matrix, acting on the state 



Pit)) = C{t) 



[2.5) 
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with Liouvillian 

m = -b 



1 - ta-nh{(3h{t)) -1 - ta-nh{^h{t)) \ 
-1 + tanh(/3/i(t)) 1 + tanh{f3h{t)) J ' 

In equilibrium, for constant h, C\p) = 0. This determines the equihbrium state 



P 



1 + tanh(/3/i) 
1 - tanh(/5/i) 



-I3ha 



[2.6] 



(2.7) 



2cosh(/5/i) 

with I 1 ) being the two dimensional unit vector. The temporal evolution of the spin is 

&,{t) = a,C{t) - C{t)(T,. (2.8) 

In addition to the spin operator a response operator az(t) is introduced. It describes 
the action of a variation 6h(t) of the external field at time t 
dC{t) 



(Jr. 



dh{t) 



[2.9) 



Acting on an equilibrium state these operators obey a fluctuation-dissipation theorem 
(FDT) 



P 



13 a, 



p). (2.10) 

For time dependent field h{t) the Liouvillians for different times do not commute 
in general. This can formally be overcome introducing a time ordering operator T 
rearranging products of objects at different time such that they are ordered from left to 
right according to decreasing time. This allows to write for t > to 



Pit) 



T 



■ r' dt'c{t') 

6-^*0 



p{to 



(2.11) 



3. Ising model on a Bethe lattice 



3.1. Dynamics 

We consider Ising spins (Tj on a Bethe lattice with couplings Jij. The Hamiltonian 
(energy) is 

H = -l^Jijaiaj -^hiai. (3.1) 

The couplings J^j are independent stochastic variables taken from some distribution 
P(J) such that Jij = and Jfj = P. The underlying vector space for the complete 
system is spanned by the direct product of the two dimensional vectors associated with 
each of the spins cTj. 

The dynamics is ruled by the Liouvillian 

C{t) = Y.C,{h) (3.2) 

i 

where Ci{ki) given by (2.6) acts on the corresponding subspace of site i and depends on 
the effective field 

dH 



(3.3) 
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Expectation values of time dependent observables with initial condition at to are written 

as 



A{t)B{t') 



N 



B{t')r 



N 



(3.4) 



■■■T 

for t > >■■■> to- 
Adopting in the following Ito-calculus the effective fields (3.3) are retarded, i.e. 

3 

It is convenient to integrate over the effective fields ki{t) in (3.4) and to take (3.5) 
into account by introducing appropriate ^-functions. Those are written in their Fourier 
representation integrating in addition over imaginary auxiliary fields ki{t). This leads 
to a path integral representation 

A{t)B{t') ■■■) = \[jv{hh} e^' 4'^ ''^'^ ^'-''^'^^ 

i 

„I]i Str,^^ {^'(*) Hi ■hj(Tj{s-)+Ci(ki{s),s)} 



X 1 r 

\n 



A{t)B{t') 



N 



(3.6) 



A state I p )^ of the system with sites is in general a superposition of direct 
products of two dimensional vectors each representing a single spin. The unit vector for 
N sites 1 1 )Ar is a direct product of single site unit vectors. ( 1 Ia? is its adjoint. 
The equilibrium state, for given J, can be expressed as 



P 



N 



N 



It obeys 



N 



0. 



(3.7) 
(3.8) 



3.2. Effective single site evolution 



For time dependent expectation values involving the spin at a single site, say site o only, 
(3.6) can be rewritten introducing an effective retarded time evolution. It is obtained 
by performing the expectation values with respect to all spins except (To- Assuming 
equilibrium initial conditions expectation values of quantities at a single site, say o, can 
be written as (3.6), 

(a„(t)---) = /p{^,A:}eV^'^^^^'-'^^^^ 



X 1 



r 



ao(t)---3^o({^}K};to)e^'o'^^°('(^)'^) 



Po{k{to))) . 



(3.9) 



The action of the remaining part of the system is contained in 

i 



x(. 


r 




7V-1 



{h-ki{s)} / ^^'.Joi J^^ds[k{s)<Ji{s )+<j{s )ki{s)] 



N-1 

(3.10) 



'N-1 
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and the initial condition in (3.9) is 

p(k)) = Ittttt e^'"'" l)- (3.11) 

'Jo 2cosh(/3A;) )o ^ > 

In the above expression the primed products and sums run over i ^ a. The expectation 

values in (3.10) refer to the system without the spin at site o. 

It should be pointed out that the above expression holds for given values of the 

couplings Jij and the average over Jjj has still to be evaluated. 

3.3. Equilibrium 

The above path integral comprises an integration at the initial time to- This is 
investigated in the following. The single time (static) expectation value of (Jo{to) is 

^o(to)) = [^S'^->'}( 1 I aUk) e^'^l 1 ) (3.12) 



2-11 

and the effective action (3.10) reduces to 

yo(^) = (e^S^^-"»)^_^. (3.13) 
In analogy to the cavity method we define Co and ko writing 

J 2tt 

This yields for the partition function and the local magnetization 

Zo = 2e'^^°cosh{pko) and {ao) = tanh{pko) (3.15) 

The two quantities defined above result in 

c„ = i{ In ( (e«E: ^..^.)^_^ ) + 1„ ( (e-^E: ^..».)^_^ ) } (3.I6) 

= A + In ( (e-'E; .'»".)^_^ ) - ,n ( (e^^El ^..".)^_^ ) }, (3.17) 

For the SK-model and for the Bethe lattice the expectation values for sites i can 
be assumed to factorize. For Ising spins and using (3.15) we find the identity 

^g±/3j„.a,\^^^^ = cosh{p Jo^) ± tanh{pki) smh{p Joi) (3.18) 
which allows to express Co and ko in terms of the fields of the adjacent sites: 

Co = ^'c{Joi,ki), ko = h + ^'u{Joi,ki), (3.19) 

i i 

c{J, k) = ^\n{ cosh.\(3J) - tanh^ipk) smh\(3J)), (3.20) 
u{J, k) = ^ artanh( tanh{pk) tanh(/3J)). (3.21) 
In particular for J ^ 0, i.e. for the SK-model, 

u( J, k) = J tanh(/?fc) c( J, k) = ^J^^2^^f^y (3.22) 

In view of this limit correlation-functions are later defined with rescaled quantities 

U{J,k) = ^ju{J,k). (3.23) 

where I is the typical size of the couplings. This ensures e.g. that the Edwards- Anderson 
order parameter stays finite in the SK-limit. 



Dynamics of an Ising Spin Glass on the Bethe Lattice 



7 



4. Slow Dynamics 

4.I. Distribution of slow fields 

Dynamics on sliort time scales takes place within a single valley of the energy landscape, 
dynamics on long time scales is supposed to be due to transitions among different valleys. 
Assume there exists a time scale t* S> 1 separating fast and slow motions. This time 
scale might be realized by some waiting time, slowly changing couplings or other means 
[4, 9]. The local field k{t) and the conjugate field k{t) are split into fast and a slow 
contributions 

k(t)^k(t) + K(t/t*), kit) ^ kit) + -kit If). (4.1) 

t* 

The fast parts are due to fluctuations around the quasi equilibrium state within a single 
valley. The slow part k{t) acts like an external field and the quasi equilibrium state 
follows this field adiabatically. Evaluating time dependent correlation functions of the 
spin (To on this long time scale, this spin can be viewed as being in equilibrium in the 
slow field resulting in 

((To(rr) (To(rr') ■ ■ ■) = tanh (/^^^(r)) tanh {(3ko{t')) (4.2) 

The effective field is given by (3.19) with k replaced by k(t), i.e. 

/€o(r) = h + Y,'u{Jo^,K.i{T~)). (4.3) 

i 

In the low temperature regime of a spin glass this field is, however, distributed and (4.2) 
has to be averaged over some distribution Po({'^o})- With (4.3) this distribution can 
be calculated from the distribution of slow fields at the neighboring sites i. Continuing 
this mapping to outer shells of the lattice an iteration scheme is set up similar to the 
course of action in the cavity approach [13]. This program involves the following steps: 
Performing the average over the couplings Jij and defining 

1 ; 

Qo{{k,K}) = \ (4.4) 

\ In-i 

this distribution is obtained by integration over k 

Vo{{k}) = |p{«:}e^-^M^'^-'^M>Q„({K,K}). (4.5) 

The average on the right hand side of (4.4) involves a corresponding joint distribution 
Qn-1 of the fields ki and Hi at the K + 1 neighboring sites. The second term in the 
exponent is the contribution of site o to the local field at site i, corresponding to the 
Onsager reaction field. 

Assuming factorization of the expectation values on the right hand side of (4.4) 
leads to 

Qoi{k, ^}) = 11' I ^{'^i' '^^} e^'' Q^{{ki, K,}) 

i 

~J O' 



^dT{k{T)u{Joi,Ki{T~)) + iii{T)u{Joi,K(T~))} ^4 
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The functional Qo{{k, k}) contains the action of all + 1 spins surrounding site o. 
Qi{{ki, Hi}), on the other hand, contains only the action of the K spins on the outgoing 
branches originating at site i. The action of the spin o"o is taken into account by the 
bond averaged exponential in (4.6). The functional Qi is given by a similar average over 
the distributions on the next shell of sites on the tree. Iterating this process, a fixed 
point distribution is assumed to exist. It is given by 



Qi{k,K}) 



-J 



X ei dr{K(r)n( J,A(r- ))+A(r)n( J^lr" ))} 



K 



(4.7) 



The functional Qo{{k, k,}) for the central spin at site a is given by the same expression 
with K replaced hj K + 1. 

The above line of arguments follows pretty much the steps in the cavity method 
calculation of Mezard and Parisi [13]. The functional fixed point equation (4.7) is, 
however, more general. The only assumptions made are the separation of time scales and 
the factorization of expectation values on different subtrees. In particular it contains 
in principle contributions corresponding to full replica symmetry breaking solutions. 
Within this formulation the replica symmetric solution is obtained by leaving out the 
reaction term in (4.7) as shown later. 



4 ■ 2. Correlation-response- functions 

There is little hope to find solutions without relying on approximations or expansions. 
In the following we investigate an expansion in powers of the difference e = -^fe^ between 
transition temperature Tc and the actual temperature T. Such an expansion is expected 
to apply in the neighborhood of a continuous freezing transition where the Edwards- 
Anderson order parameter qea ~ ^ near T^. This is the case for the SK-model. 

Rather than working with the full functional Q{{k, k}), it is sufficient to investigate 
its moments, the correlation-response-functions on the long time scale. With the rescaled 
quantities U{J,k), (3.23), they are defined as 

Cnmin, ■ ■ ■ , r„; r(, ■ ■ ■ , r^; J) =1 V{k, k} e/d'^^(-)i'^-''W} Q{{k, k}) 

n m 

y<X{U{J,K{r,))X{k{rl). (4.8) 

Inserting (4.7) and expanding the exponential, a hierarchy of equations for the 
correlation-response-functions is obtained which is derived along the following lines: 
Evaluating a function Cnmiji) ■ ■ ■ , T"n; t{ , ■ ■ ■ , r^; J) there is the product over U (J, i^{Ty)) 
involving k at the times r^. The factors k(t^) can be replaced under the path integral by 
functional derivatives 5/5k{t^) acting on Q{{k, n}), at least as long as all time arguments 
are different (otherwise the derivatives could act on the products of f/(J, k(tj^)) as well). 
The derivatives acting on Q create additional terms involving k(t^). Expanding the 
second term in the exponential of (4.7), internal integrations over times dp are generated. 
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They contribute additional factors now depending on /t(o"p). In order / of this expansion 
there are altogether functions of + m + / discrete times. For r ^ r^, r^, the path 
integral over k{t) can be performed resulting in ^(r) = 0. This leaves integrations over 
k and k at those discrete times only. With (4.7) they reduce to functions of correlation- 
response-functions. This general strategy is applied in the following to correlation- 
response-functions of lowest order. 



Tc-T 



For the SK-model C„ 



^ g(n+3m)/2 



5. Expansion around Tc 

5.1. Leading order for general P{J) 

We investigate an expansion in powers of e 
is found [3, 4] and the same scaling with e is proposed for the present investigation as 
well. Later it is shown that this scaling is fulfilled in a consistent manner. For simplicity 
the following discussion will be restricted to h = 0. 
We may start with an expansion of (4.7) 

Q{{k, /€}) ^ 1 + dadcr' k{a)k{cr') 6*2,0 ((X, (r';J'f 

+I'^J dada' k{a) U{J' , K{a')) Ci,i{a, a'-J'f 

ji 

d(Tidcr2do-3dcr4 k{ai)k{a2)k{a-i)k{ai) C4,o(o-i, c^2, 0-3, cfa] J') 

ji 

daido-2d(T3d(T4 k{ai)k{a2)k{a:i) U{J\ «:((T4)) C3,i(cri, 0-2, crz] o^; J') 

+ ■■] ■ (5.1) 
Evaluating 6*2,0 (t, t'; J) and keeping terms ~ 
C*2,o(r,r'; J) = «:} e^'^^('^H'^-''('^»t/( J, «:(r))f/( J, ^(r')) 



^I^Kldada' k{a)k{a') C2,o{a, a'; J')" 



+^I'^kJ dcrid(J2d(J3dcr4 k{ai)k{a2)k{as)k{a4) 



-J' 



J'- 



(5.2) 

rJ" 



C*4,o(cri,cr2,o-3,o-4;-^') +3(A'- l)C*2,o(o-i,a2; J') 6*2,0(0-3, (74; J" 

Performing the path integration over K{a) for a t,t' the integrals over a, reduce to 
sums over a = r, r' and the path integrals become ordinary integrations over variables 
k = k{T), K = k,{t), k' = k{T') and n' = k(t'). 
It is convenient to introduce vertices 

K(J) = ^f/(J,«:)L=o (5-3) 



with = for odd n. In the following we use the notation A{J)B{J) AB. Then 
(5.2) yields 

C*2,o(r, r; J) = \l^KV^{J) + \l^K{K - l)V2{J)V^{J) 



C2,o(r,r) + C*2,o(r',r') 
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y-C^o{r. r') + ^I^K V2{J)V^{J) [ CV^(r, r, r, r') + CV^(r, r', r', r') ] . (5.4) 
Evaluation of C/^fl, again with /i = 0, in order results in 

C4,o(ri, r2, rs, J) = l'KV,\J) C^(ri, ra, rg, + l'K{K - 1)V,\J) (5.5) 

) + C2 ,o(ti, T4) C2 ,0(^2, T3) 

and performing the average over J we get 



C4,o(n,r2,r3,r4)=/^ir(i^-l) 



C2,o(T17 ''"2) C'2,o(''"3) T4) 



+ (^2,0(^-1,73) C2,o(t2,T4) + C2,o(Tl,r4) 6*2,0(7-2, T3) 

With this C2,o becomes 
C2,o(r,r';J) = 



(5.6) 



C2,o(r',r') 



xC2,o(r,r'). 



(5.7) 



It should be noted that the correlation function C4 o has been broken down completely 
to products of pair functions 6*2,0- This can be done for correlation- response-functions 
Cn,m -with general n and m in higher orders as well. In this case the products contain 
also response functions Ci 1. 



Averaging over J a non trivial solution with C2,o(t, t') 7^ requires that the 
expression in the curly bracket on the right hand side equals 1. Including terms of higher 
orders in e shows that a corresponding expression holds only for the time derivative. 
Assuming time translational invariance we can write 

(5.8) 

with 



drC2,o{r)=Y{T)drC2,o{r) 
Y{t) = I^KV^ + I^K{K - 1] 



C2,o(0). 



(5.9) 



In first order Y{t) is constant and (5.8) would hold for C2,o(t) as well. Including higher 
orders, however, it depends explicitly on r and (5.8) holds for the derivative only. A 
non trivial solution requires Y{t) = 1 for all r and this eventually determines the full 
time dependence of 6*2,0 (t). 

The evaluation of C*i.i follows similar lines. Since C*i,i ~ e 
kept, and the response function becomes 

C*i,i(r, r'; J) = I'KV2{J) ^^(r, r') + '^l'K{K - 1)V,{J) C^(r, r) T^(r, r' 



^ contributions ~ are 



1 r4 



I'K{K - 1)V2{J) V4C*i,i(r; r') C*2,o(r',r') 
+|,/^K\/4(J) W^(r,r,r;r') + II^KV2{J)W^,{t,t',t';t'). (5.10) 
Likewise the four point function on the right hand side obeys 
C*3,i(ri, r2, T3; n; J) = I^KVi{J) m7(ri, ts, t,- t,) 



+l''K{K-l)Vi{J) Wi,i(n;r4) C*2,o(r2,r3) 



+ V2Ci^i{T2]n) C*2,o(ri,r3) + V2Ci^i{T^]n) C*2,o(ri,r2) 



(5.11) 
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The average over J requires some attention because on the right hand side of (5.10) and 
(5.11) the averages contain different vertices Vn{J)- Nevertheless they can be evaluated 
successively. Multiplying (5.11) with V2(J) and performing the J-average, ¥203^1 can 
be expressed in terms of V2Ci,i and V2Ci,i. Multiplying (5.11) with V4(J), averaging 
and using the result for V2C3,i, ^46*3^1 is again written in terms of pair functions. This 
is inserted into (5.10) resulting in 



Ci,i(r, r , J) = ^/ KV2{J) + K[K-l) 1 - I^KV^ 

xT^(r; r') + \I^K{K - 1) V2( J)"C^(r , t') T^(r; r'). (5.12) 
This yields in lowest order 



Wi,i(r; r') = PK V^V^ V2Ci,i(r; r'). (5.13) 



Inserted into (5.12) a closed equation for V2C\^\ is obtained 



WM(r;r') = [I'KV:^ + \1'K{K - 1) ^ ^f,"—, C2,o{r, r) (5.14) 



+\I^K{K - 1) 
which is again of the form 



C2,o(r,r')}Wi,i(r;r') 



Wi,i(r) = F'(r)V2Ci,i(r) (5.15) 

assuming time translational invariance. 

Pointing out the difference between Y{t), (5.9), and Y'{t) we write 



Y\r) = F(r) + l/^i^(K-l){ ''^^^^ ^^^^^ - [l-J^i^W] 1^)^(0) • (5.16) 



This allows for the following non trivial solutions: 



(a) F(r) = 1, Y'ir) ^ 1 : C2,o(r) ^ 0, Ci,i(r) = 0. 



(b) F(r) ^ 1, Y\r) = 1 : C2,o(r) = 0, Ci,i(r) ^ 0. 



(c) Y{t) = Y'{t) = 1: C2,o(r)^0, Ci,i(r) ^ 0. 

A solution of type (a) has been found for instance for the spherical version of the 
SK-model[16]. It does not describe glassy behavior since the response function vanishes 
on the long time scale, although the correlation function shows non trivial properties. 
A solution of type (b) does actually not exist as shown later. 

Only the last solution is supposed to be characteristic for a spin glass phase. Since 
PKV2 = l + 0(e) the last term in (5.16) does not contribute in order e and only the first 
term in the bracket has to be taken into account. This means, however, that solution 
(c) requires 



Vi ViVi = ViXWi. (5.17) 
With (3.21), (3.23) and (5.3) 

V2{J) = tanh(/3J) Vi{J) = -2/52J-1 tanh(/3J) (1 - tanh'^{/3J)) (5.18) 
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and 

viij) viij')v,ij') - vi{j) V2{j')v^{j') 

= 2(3^rHanh\(3J)tanh\/3J'){tanh\/3J) -t8.nh\/3J')f. (5.19) 

This shows that (5.17) can be fulfilled only with J = ±1 and that this kind of solution 
does not exist for more general P{J). This result is quite remarkable, indicating the 
special role of the binary distribution of couplings. 



5.2. Next to leading order for J = ±7 

The correlation-response-functions (4.8) have the symmetry 

Cnmi- ■■;</) = i-irCnmi- " • ; "J). (5.20) 

For binary couplings J = ±/ the bond average therefore becomes almost trivial. For 
simplicity we set h = in the following and assume time translational invariance on the 
long time scale. The strategy to evaluate contributions of higher orders in e is analogous 
to what has been described in section 5.1. Evaluating 6*2,0 in order requires to compute 
6*4,0 and C*6,o in the same order and to break them down successively to products of pair 
functions C*2.o- In addition in this order an internal integration shows up taking into 
account the third line on the right hand side of (5.1) leading to a contribution of the 
form 

C*2,o(r, T') = --- + V2jda C*3,i(r, r', a; a+). (5.21) 

The four point function C*3,i can again be broken down to products of pair functions 
according to (5.11). 

Using the notation 

q(r - t') = C*2,o(r, r'; J), f(r - r') = C*i,i(r; r'; I), K = V^) 

and collecting all contributions up to order one obtains 

V2V4 



(5.22) 



g(r) = I'K g(r) + rK{K - 1) 



+I^K{K - 1) 



1 



1 - I^K 



X 



\{V2Ve+V^)+3I^K 



-2 + 3I^K{K 



^2 



V2 V4 



Vf + 2rK 



V2 V4 



-I^K{K 



d(T|r(r -|- cr)g((T) -|- g(r — cr)r((T)|. 



1 - I^KVi 

The corresponding calculation for the response function yields 



r{r)\-q{r) 

(5.23) 



f(r) = I'K f{T) + rK{K - 1] 



g(0)f(r) 



+I^K(K - 1) \ ^{k -2 + 3I^K(K - 1) ^ j\ 
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HV2Ve + Vi)+3I^K 



2 ''4 



I r\n- <r( ' 

1 - I^K 



2 , or4r^ ^2 



+I^K{K-1)- Jl,^ /daf(r-a)f(a). (5.24) 



5.3. Ultrametric time parametrization 

The low temperature phase of the SK-model is characterized by a hierarchy of long time 

scales ranging from t* to some longest time scale t^. Eventually the limit — > oo is 
taken. In order to keep track of the long time scales it is convenient to introduce the 
parameterization [4] 

t = tl-^^\ x{t) = l-^ (5.25) 

In(t^) 

and write 

g(r) = Q(x(t*r)), f(r) = -t* x{t*T) R{x{t*T)). (5.26) 
With the above definition 

x{t + = 1- , + ^ = x{t) - ^"^'7- . > (5.27) 

In(t^) In(t^) 

and for t > t' and with it x(t) < x{t') 

InM + ^^(*)-^'{*')^ fx{t)-x{t') 

x{t + t') = x{t) - ^ I - ^ ^ x{t) - \—- ^ xit). (5.28) 

In(t^) In(t^) 

The corresponding result for t < t is x{t+t') ^ x{t'). This yields the ultrametric relation 

x{t + t') = x{ max(t, t')). (5.29) 
With it the integrals in (5.23) and (5.24) become 
J da f{T — a)q{(j) + J da ^f{T + a)q{a) + q{T + a)f{a)^ 

2Q{x) Cdx R{x) + 2 / dx Q{x)R{x') (5.30) 

Jx Jo 

and 

rda f(r - a)f{a) -2t*x{t * r) R{x) / W i?(a;') (5.31) 

Jo Jx 

This allows to write (5.24) in the form 

R{x) =Y{x)R{x) (5.32) 

with 

Y{x) = Yo + YiQ{l) + ^^20^(1) + iy2xQ'(x) + Yr /W i?(x') (5.33) 

Jx 

and the coefficients 

yo = /^i^^^ (5-34) 
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- - 1) r^^{K - 2 + si^m - d^^^} 



^{Vi + 2I^K ^ (5-37) 

and 



1 - f ' 



= 2I'K{K - 1) (5.38) 
The correlation function (5.23) becomes 
Q{x) = {Fo + Y^Q{1) + \Y2Q\l) + ;|F2,g'(x)|Q(a:) 

+n|g(x)j^dx'i?(x') + j\x' Q{x')R{x')^. (5.39) 
Differentiation with respect to x results in 

a,g(x) = F(x)9,Q(x) (5.40) 
with the same Y[x) as above. This allows a non trivial solution for Y[x) = 1. 

5.4- Results in next to leading order 

In particular for x = 1 and in second order we get 

Lowering the temperature a non zero Edwards- Anderson order parameter qea > shows 
up first at a critical temperature where 

Yo = I^K Vi = K tanh\(3I) = 1. (5.42) 

This determines the critical temperature 

Tc = ^ (5.43) 

artanh(l//^) ^ ^ 

and choosing 

/ = artanh(l//^) (5.44) 

the critical temperature is Tc = 1. 

Expanding the vertices and coefficients (5.36) to (5.38) in powers of e = 1 — T 

Fo = 1 + + I {IK -31 + 2VK)^e\ 

Yi = -2I^{K - 1) - AI^{iVk + 1){K - l)e, 

Y2 = lOI^iK'^ - 1), Fax = 4:I^iK^ - 1), Yr = 2IVk (5.45) 
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and (5.41) becomes 



1 

= p + 2- , 

K 

Differentiating (5.33) with respect to x gives 



e\ (5.46) 



ir y K 



This verifies that R{x) ~ e^. 

In equihbrium a fiuctuation- dissipation theorem (FDT) holds, i.e. r{t) = —f3dtq{t). 
As pointed out earher this is also expected to hold in the glassy state for t < t*. On the 
long time scale, however, the FDT is violated. We may introduce a measure X for the 
violation of the FDT 

f(r) = -f3 X{t) drqir), R{x) = f3X{Q{x)) 9,g(x) (5.48) 

with 

Xiq) = ^q = 2I^^^^^q for < q < qEA- (5.49) 

Ir V K 



in lowest order. 

The fact that X depends on Q only holds in higher orders as well. For t < t* the 
correlation function obeys q(t) > qEA and we can set X{q) = 1 for t < t*. This quantity 
is related to Parisi's overlap distribution function P{q) = dqX{q). 

Eq.(5.47) does not fix the actual form of Q{x) or R{x). Using the so called Parisi 
gauge [17], l3X{Q{x)) = x, 

Yr 

Q{x) = TF~x for X < X* (5.50) 



with 



x* = ^qEA (5.51) 

1 



and t* = . 

For K ^ 00 the known results for the SK-model are recovered. 



6. Stability analysis 

6.1. Solutions for general distribution of couplings 

For general distributions P{J) the difference Y'{x) — Y[x), given in (5.16) is non zero. 
In order e the last term ~ 1 —PKV2 does not contribute. The remaining term contains 

TfWv^-WW^>o, (6.1) 

which is positive according to (5.19) unless J = ±1. Accordingly Y'{x) > Y{x) for 
general distributions P{J). 

This still allows for solution (a) of section 5.1 with R{x) = 0. The last term 
in (5.33) vanishes and Q{x) = qEA, given in (5.46). This solution does not show any 
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time dependence on long time scales. This property is shared by correlation functions 
of higher order, C„^o(ti ■ ■ ■ Tn) = Cnfl- Correlation-response-functions Cn,m{- ■ ■) with 
m > vanish. Reconstructing the full distribution Q{{k,K,}) from (5.1) one finds the 
replica symmetric solution [13] 



Q{{k, k}) = jdRP{R) e/d-«(-)W-)-^l 



with 



K 

P{k) = \{jdKi P{Ki) 5{h + ^ n( Ji, 

j=l i 



Ki ] - K] 



In next to leading order 

Y{x) = Yo + Y, qea + + |rF2x] qIa 

which yields 



Qea 



13 + 1 



3 K I,/k\ 
As shown below this solution is unstable. 

Solution (b) of section 5.1 requires Q{x) = and with (5.33) 



d^Y{x) = Yo + Yr dx' R{x') = 1. 

Jx 

This has, however, no solution with R{x) 7^ 0. 



(6.2) 
(6.3) 
(6.4) 
(6.5) 



(6.6) 



6.2. Stability analysis 

Various criteria can be used to test the stability of a given solution. One may for 
instance ask whether the correlation-response- functions on the short time scale approach 
the values determined by the dynamics on the long time scale. For example the pair 
correlation function q{t) = {a{t)a{0)) should approach qea for t ^ t*. This decay is 
expected to be algebraic, but a corresponding analysis is outside the scope of this paper. 

The general procedure used in statics as well as in the present formulation of 
dynamics consists in iteratively connecting subtrees to a new tree and expressing the 
properties at the new vertex by properties at the base vertices of the subtrees searching 
for fixed points of this mapping. A necessary condition for stability is the decay of small 
perturbations under this mapping. Eqs.(5.23) and (5.24) are actually such mappings. 

For T > Tc the trivial solution Q{x) = and R{x) = is expected to be valid. A 
small perturbation 6Qi{x) on site i is mapped onto 

SQo{x) = PKV^5Qi{x) (6.7) 

at site o, and a corresponding mapping for Rq{x). Since PKV2 < 1 for T > Tc the 
trivial solution fulfills this stability criterion. By the same token this solution is unstable 
for T < Tc. 

The stability criterion for solution (a) reads with (6.4) and Y{x) = 1 

6Qo{x) = [1 + iFsx qIa] SQ^i^) > SQii^) (6.8) 
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and the same for 6R{x). This solution is therefore unstable for T < as well. 

Finally solution (c) for J = ±/ is investigated. We may test the stability with 
respect to a perturbation 6Q{x) and SR{x) = X[Q[x)) dxSQ{x). This yields with 
Y{x) = 1, 6Qo{x) = SQi{x) and dRo{x) = dRi(x). This means that this solution 
is marginal with respect to perturbations of this kind. Investigating a perturbation 
6R{x) with 6Q{x) = one obtains 



Considering the second contribution a perturbation at some value x' creates fluctuations 
at X < x' only. This means that the perturbation at x' is not enhanced due to this 
term and this solution is marginal with respect to a perturbation of this kind as well. 
Such a marginal stability criterion is actually expected because of the reparametrization 
invariance mentioned at the end of Section 5.4. 

7. Discussion 

In this paper we have shown that the long time dynamics of an Ising spin glass with 
binary couplings J = ±/ on a Bethe lattice is of the form known from the Sherrington- 
Kirkpatrick model. This is remarkable insofar as the interaction is restricted to nearest 
neighbors. In contrast to a lattice in finite dimensions, however, the typical size of closed 
loops scales with In(iV) for a system of sites. 

The second main result of this paper, the failure of a corresponding solution for 
general distributions P{J), is unexpected, and the nature of the low temperature phase 
in this case is not known. A similar breakdown of a replica or cavity method calculation 
for general P{J) can not be excluded since the relevant equations have been evaluated 
for binary coupling only [13]. It might be of interest to perform an expansion of the 
relevant equations around Tc, possibly within the extended replica scheme proposed by 
de Dominicis et.al [17]. In this scheme two order parameter functions corresponding to 
Q{x) and R{x) are used. 

This paper is dedicated to David Sherrington on the occasion of his 65th birthday. 
References 

[1] M. Mezard, G. Parisi and M.A. Virasoro, Spin Glass Theory and Beyond (World Scientific, 
Singapore, 1987) 

[2] G. Semerjian, L. F. Cugliandolo and A. Montanari Journal of Statistical Physics, 115, 493, (2004) 
[3] H. Sompolinsky and A. Zippelius,Phys. Rev. B 25, 6860 (1982) 
[4] H. Horner, Z. Phys. B 57, 29, 39 (1984) 

[5] D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 1972 (1975) 

[6] T. R. Kirkpatrick and D. Thirumalai, Phys. Rev. Lett. 58, 2091 (1987); Phys. Rev. B 36, 5388 
(1987) 

[7] A. Crisanti, H. Horner and H.-J. Sommers, Z. Phys. B 92, 257 (1993) 
[8] L. F. CugUandolo and J. Kurchan, Phys. Rev. Lett. 71, 173 (1993) 




(6.9) 



Dynamics of an Ising Spin Glass on the Bethe Lattice 



18 



[9] M. Freixa-Pascual and H. Horner, Z. Phys. B 80, 95 (1990) 

[10] A. Barrat, R. Burioni and M. Mezard, J. Phys. A 29, L81 (1996) 

[11] H. Horner, Z. Phys. B 86, 291 (1992) 

[12] H. Horner, arXiv:0707.2714, to be pubhshcd in Eur. Phys. J. B (2008) 

[13] M. Mezard and G. Parisi, Eur. Phys. J. B 20, 217, (2001) 

[14] P. Mottishaw, Europhys. Lett. 4, 333 (1987) 

[15] C. De Dominicis, Y.Y. Goldschmidt, J. Phys. A 22, L775 (1989), Phys. Rev. B 41, 2184 (1990) 

[16] W. Zippold, R. Kahn and H. Horner, Eur.Phys.J. B 13, 531 (2000) 

[17] C. Dc Dominicis,M. Gabay and H. Orland, J. Physique Lettres 42, L523, (1981) 



